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I.   INTRODUCTION 

During  the  past  few  years,  there  have  been  many  advance- 
ments in  numerical  techniques  for  predicting  the  behavior  of 
fluid  flows.   For  example,  several  computer  models  have  been 
developed  by  Gosman,  Spalding  and  others  [1,2,3]  which  use 
the  mass,  momentum  and  energy  conservation  equations  reduced 
to  finite  difference,  nonlinear  algebraic  form.   The  develop- 
ment of  reliable  computer  programs  of  this  type  greatly  bene- 
fits engineering  analysis  in  such  widely  varying  fields  as 
meteorology,  aerodynamics  and  gasdynamics. 

The  earlier  two-dimensional  computer  codes  were  based  on 
vorticity  (qj)  and  stream  function  (<J/)  [1,2,5].   This  form  of 
the  governing  equations  eliminates  pressure  and  velocity  from 
immediate  consideration.   Pressure  is  normally  calculated 
only  after  a  converged  solution  is  obtained.   This  technique 
has  several  inherent  disadvantages: 

1.  It  results  in  large  errors  in  the  predicted  pressure 
distributions  in  all  but  quiescent  flow  regions  due  to  the 
higher  order  dependence  of  the  pressure  gradient  on  stream 
function  [6 ] . 

2.  It  is  usually  restricted  to  constant  density  flows 
or  to  flows  in  which  density  varies  only  with  temperature 
[3,6]. 

3.  The  boundary  conditions  are  difficult  to  specify 
[3,5] 

4.  Considerable  difficulty  has  been  experienced  in 
arriving  at  converged  solutions,  especially  for  non-uniformly 


spaced  grids  and  high  flow  rates  [2,5,6,7]. 

5.   The   ip— oj   model  is  not  easily  extended  to  three- 
dimensional  flows  [3], 

To  overcome  these  difficulties,  emphasis  has  been  placed  on 
developing  computer  codes  based  on  velocity  and  pressure,  the 
primitive  variables. 

A  major  problem  with  any  new  computer  model  is  model 
validation.   The  difficulties  of  collecting  accurate  empirical 
data  are  multiplied  when  investigating  three-dimensional  and/ 
or  reacting  flows.   In  addition,  many  variables  within  these 
flows  are  not  readily  measurable  (turbulence  intensities,  etc.). 

An  effort  to  utlize  elliptic  computer  models  which  can 
handle  turbulent,  reacting,  variable  density  flows  at  high 
subsonic  and  sonic  velocities  has  been  underway  at  the  Naval 
Postgraduate  School  for  several  years.   Two  specific  areas 
which  have  been  investigated  are  flows  in  a  turbojet  test  cell 
and  in  the  combustion  environment  of  a  solid  fuel  ramjet. 

It  is  important  to  have  the  capability  to  test  high  per- 
formance jet  engines  throughout  their  operating  envelope 
under  conditions  which  approximate  installed  conditions.   This 
is  often  accomplished  in  blockhouse  type  installations  called 
turbojet  test  cells  (TJTC) .   The  typical  test  cell  incorporates 
an  inlet,  a  horizontal  test  section  and  a  vertical  exhaust 
stack.   The  engine  to  be  tested  is  normally  mounted  near  the 
center  of  the  cell  to  allow  the  development  of  a  nearly  uniform 
engine  inlet  velocity  profile.   The  engine  exhausts  into  an 
augmentor  tube  which  entrains  additional  air  for  exhaust  gas 
cooling  and  dilution.   The  quantity  of  this  secondary  air  is 


crucial  to  proper  engine  testing  and  test  cell  performance. 

Testing  today's  high  power  and  high  mass  flow  engines  in 
these  installations  produces  a  myriad  of  noise  and  air  pollu- 
tion problems.   Cell  modifications  must  often  be  made  to  mini- 
mize these  problems.   This  fact  coupled  with  the  future  need 
for  larger,  more  expensive  test  cells  to  replace  obsolete  cells 
and  to  accommodate  new  generations  of  high  technology  engines, 
makes  the  development  of  reliable  modeling  methods  imperative. 
The  frequently  used  one-dimensional  models  are  not  adequate 
for  predicting  the  details  of  the  complicated  flows  within  a 
turbojet  test  cell  and,  therefore,  the  cells  often  do  not  per- 
form to  their  designed  limits.   An  accurate  flow  model  would 
provide  a  needed  design  tool  which  could  help  prevent  costly 
design  errors  and  improve  operating  efficiency. 

A  two-dimensional   ^-co   computer  code  was  used  to  analyze 
the  flows  in  a  full  scale  and  a  subscale  turbojet  test  cell 
at  the  Naval  Postgraduate  School  [2,6].   Experimental  data 
from  the  subscale  test  cell  have  been  compared  with  computa- 
tions made  with  this  computer  model  [6].   A  primitive  variable 
(u-v-p)  computer  model  could  improve  this  capability  by  ex- 
tending it  to  specific  geometries  and  flow  rates  that  the 
tJ>-w   model  was  incapable  of  predicting.   In  addition,  a  u-v-p 
model  would  more  readily  allow  variable  density  flows  to  be 
analyzed  and  should  more  accurately  predict  augmentor  pressure 
distribution.   An  application  of  a  primitive  variable  model 
to  jet  pumps  has  been  reported  by  Croft  and  Lilley  [4]. 


The  primary  objective  of  this  investigation  was  to  adapt 
a  primitive  variable  model  to  the  turbojet  test  cell  geometry 
and  flow  conditions  and  to  validate  that  model  with  data  from 
a  subscale  turbojet  test  cell.   Possible  utilization  of  the 
model  for  predicting  optimum  placement  of  sound  suppressors/ 
jet  break-up  devices  within  the  augmentor  tube  was  also  inves- 
tigated. 


II.       MODEL    OVERVIEW 

A.  INTRODUCTION 

The  computer  model  used  in  this  study  was  adapted  from 
the  CHAMPION  2/E/FIX  computer  program  developed  by  Pun  and 
Spalding  [8],   CHAMPION  is  a  TWO- dimensional  Elliptic,  FIXed 
grid  computer  program  which  provides  a  solution  of  the  conser- 
vation equations  for  recirculating  flows  in  finite  difference 
form. 

B.  ASSUMPTIONS 

The    flow  was    assumed   to  be   steady,    two-dimensional    and 
subsonic.      For   simplicity   the   value   of   specific  heat    (C   )    was 
assumed  to  be    constant   although   its    dependence   on    temperature 
and/or   composition    could  easily  be    included. 

A  modified  Jones-Launder    [8,10,11,12]    two   parameter   tur- 
bulence  model  was    incorporated   to   calculate    the  effective 
viscosity.      It   uses    five   empirical    constants    (Table    I)    and 
requires    that  two    additional    variables,    turbulence   kinetic 
energy    (K)    and   turbulence    dissipation   rate    (e) ,    be   evaluated. 
Effective    viscosity  was    calculated   using   the    formulas: 


yeff  =   ulam  +   ^t  (1) 


where 


-> 


Ut   =   CD   p    KVe  (2) 


The    turbulent  Prandtl    and   Schmidt  numbers   were    taken 
equal    to   unity   and,    therefore,    the    turbulent   Lewis   number  was 


C.  C~  C_  a,        ~,  a        a 

1  2  D  k,eff  e,eff 


1.43         1.92         0.09  1.0  1.3 


TABLE    I.       K-£    TURBULENCE    MODEL   EMPIRICAL    CONSTANTS 


unity    [9  J.      The    laminar  Prandtl   number  was   also    taken    to  be 

unity. 

C.       GOVERNING   EQUATIONS 

The  conservation  equations  for  axi-symmetrical  flows 
with  no  tangential  variations  can  be  put  into  the  general 
form    [ 8 ] . 


!_  <pu,)  + 1 1_  (prvW  .  |_  ff#  || .  1 1_  (rr^  ||  =  s$  o) 

convection  terms        diffusion  terms     source 

terms 


where      $      stands    for  the   dependent  variable    (u,v,k,e,h,    etc.) 
being   considered    (<£   =   1    for   the    continuity   equation)  ,       r        is 
the   appropriate   effective   exchange    coefficient    for   turbulent 
flow   and      S        is    the    "source    term"     (Table    II) .      The   energy 
equation   in   terms    of   stagnation   enthalpy   has    no    source    terms 
since    the    turbulent  Prandtl    and  Schmidt   numbers   were    chosen 
as    unity    and   radiative    transport  was   neglected    [1,3].      Thus 
the    stagnation   enthalpy    is    given  by: 


h   =   h   +    (u2   +   v2)/2    +   K  (4) 


where    for  non-reacting    flows: 

h   =   C      T  (5) 

P 

The  calculation  of  temperature  was  made  using  equations 
(4)  and  (5) .   Density  was  calculated  from  the  perfect  gas  law 

p  =  P/RT  (6) 

where  r  =  r/m  (7) 
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D.  CONSERVATION   OF   MASS 

On   each    radial    line    the  mass    flow   rate  was    calculated 
using   the    local   density.      The  error  in  mass    flow    (compared 
to   the  summation  of   "mass-in"    at  all   upstream  boundaries)    was 
used  to   uniformly   adjust  the  axial   velocity  over  the   entire 
line.      This  process   ensured  that  overall   continuity  was   satis- 
fied on   the   line.      The   pressure   at  all   downstream  locations 
was    then   adjusted   to   approximately   correct   for   the   momentum 
imbalance    created  by    the   uniform  axial   velocity   adjustment. 
A   "pressure    correction"    equation  was    then   solved   for  each    cell 
on   the    line.      Local   cell    velocity    (axial    and   radial)    and  pres- 
sure were    then   adjusted   to   satisfy   cell-wise    continuity.      The 
details   of   this   procedure   are   presented  in   reference    8. 

E.  BOUNDARY  CONDITIONS 
1.   Introduction 

Fixed  boundary  conditions  were  set  at  the  desired  or 
experimentally  determined  value  and  held  constant.   Specified 
gradient  boundary  conditions  were  handled  by  setting  the 
appropriate  convection/diffusion  coefficient  to  zero  in  the 
finite  difference  equation  ("breaking  the  link")  and  then 
entering  the  appropriate  gradient  through  linearized  "false" 
source  terms  [8].   The  geometry  and  appropriate  boundary  con- 
ditions for  the  test  cell  are  summarized  in  figure  1. 
2.   Inlet 

Although  not.  a  computer  program  limitation,  "plug 
flow"  was  assumed  at  the  nozzle  exit  and  cell  inlet  planes. 
The  (secondary)  flow  inlet  of  the  TJTC  was  recessed  approxi- 
mately 0.3  meters  (figure  1)  to  allow  a  velocity  profile  to 


develop  over  the  length  of  the  engine.   Turbulence  kinetic 
energy  was  selected  to  be  uniform  with  a  value  which  corres- 
ponded to  the  approximate  turbulence  intensity  of  the  inlet 
flow. 

3.  Axis  of  Symmetry  and  Exit  Plane 

Radial  and  axial  gradients  were  set  equal  to  zero  on 
the  center  line  and  exit  respectively.  The  radial  flow  velo- 
city was  equated  to  zero. 

4 .  Solid  Boundaries 

All  solid  boundaries  were  considered  adiabatic  with 
both  velocity  components  equal  to  zero  ("no  slip"  condition) . 

For  simplicity,  a  two  part  boundary  layer  was  used. 
The  border  between  the  laminar  sublayer  and  the  turbulent 
layer  was  taken  at  y   =  11.5  [8],   y    was  evaluated  at 
each  near  wall  node  (p) , 


^+  ■  (P«/B1„)(V)1/2  (8) 


where,  for  y    >  11.5 
P  ~ 


\   -   CD1/2  P  Kp  (9) 


x   was  assumed  uniform  from  the  wall  to  the  near  wall  grid 

w  3 


point.   Thus, 


y  +  =  C  1/4  p  K  1/2  6/u,  (10) 

^p      D       p       Mlam 


If  y    >  11.5  ,  the  wall  shear  stress   (t  )   was  calculated 
p  —  w 

using  the  formula: 

10 


T      =    C      '        p    K      =    p    C      '       K/        (u/u    ) 

w  D  pD  p  ' 


=  *  S1/4  p  %  K^/m^sc^V72/^) 


(11) 


where 


u+  s  i  ln(E  y  *)  (12) 

K  p 

Wall  shear  stress  was  evaluated  for  y    <  11.5   from  the 


formula: 


T  =  w-    u  /6  (13) 

w    ^lam  p' 


Due  to  the  steep  gradients  of  properties  in  turbulent 
flows  near  solid  boundaries,  the  source  terms  for  K   and   e 
at  near  wall  nodes  were  expressed  in  terms  of  the  wall  shear 
stress  [1,8].   t    also  provides  the  boundary  condition  for 
the   u   and  v  equations.   In  the  following  equation  for  tur- 
bulence dissipation  rate  (e)  at  a  near  wall  node  (p) ,  the 
length  scale  is  presumed  proportional  to  the  distance  from  the 
wall  (6)  . 


£   =  Cn3/4  K  3/2/k6  =  K  3/2/2.436  (14) 


Equation  (13)  implies  that  the  wall  shear  stress  is  calculated 
assuming  a  linear  velocity  profile  when   y    <  11.5.   A  near- 
wall  grid  point,  therefore,  can  lie  within  the  laminar  sublayer, 

but  the  source  terms  for  K   and  z      imply  that   u  ^.c/u,  _   is 

c   -1  err   lam 


11 


much  greater  than  one  [10,11].   This  fact  precludes   y    from 
being  significantly  less  than  11.5. 
F.   SOLUTION  PROCEDURE 

Five  variables  (u,v,K,e   and  h)  were  solved  using  equation 
(3)  in  finite  difference  form.   The  line  by  line  iterative 
procedure  employed  upwind  differencing  and  under  relaxation  to 
promote  convergence  [8].   Pressure  (relative  to  a  selectable 
position  and  magnitude  within  the  grid)  was  obtained  from  the 
mass  conservation  imposed  on  each  radial  grid  line  and  on  each 
nodal  control  volume  as  discussed  above.   Effective  viscosity, 
temperature  and  density  were  also  obtained  as  described  above. 
A  more  detailed  explanation  of  this  procedure  can  be  found  in 
reference  8. 
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III.   TURBOJET  TEST  CELL  NOISE  CONSIDERATIONS 

A.   INTRODUCTION 

Lighthill  [13]  and  others  have  built  a  general  physical 
theory  of  aerodynamic  sound  production.   He  has  also  special- 
ized that  theory  to  the  study  of  turbulence  as  a  source  of 
sound  [14].   Experiments  have  shown  [14]  that  aerodynamic  noise 
generated  by  a  subsonic  cold  jet  in  the  absence  of  any  fluc- 
tuating forces  between  the  fluid  and  solid  boundaries  exhibits 
the  following  characteristics: 

(1)  The  acoustic  power  output  varies  as  a  high  power 
(near  the  eighth)  of  the  jet  velocity. 

(2)  A  very  broad  spectrum  with  a  peak  frequency  near 

2 
U  /2d  ,  where  U   is  the  jet  velocity  and  d   the  nozzle  dia- 
meter. 

(3)  Almost  all  the  sound  is  radiated  in  directions  making 
an  acute  angle  with  the  jet  axis. 

(a)  Higher  frequencies,  which  are  apparently  emitted 
mainly  from  the  vicinity  of  the  nozzle  exhaust  plane,  are  ra- 
diated at  an  angle  of  45  ,  or  slightly  less,  to  the  jet  axis. 

(b)  Lower  frequencies,  which  are  apparently  emitted 
from  the  jet  five  to  twenty  diameters  downstream  of  the  nozzle, 
are  radiated  at  much  smaller  angles,. 

Lighthill  postulated  [14]  that  the  high  frequency  sound 
emanates  from  the  heavily  sheared  mixing  region  near  the  nozzle 
exhaust  and  that  the  low  frequency  emanates  from  the  region  of 
more  nearly  isotropic  turbulence  farther  downstream  in  the  core 
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of  the  jet. 

While  a  myriad  of  research  has  been  done  on  the  generation 
and  propagation  of  noise  from  a  free  jet,  little  research  has 
been  directed  toward  the  problem  of  aerodynamic  noise  produced 
by  confined  jets.   The  propagation  of  sound  associated  with 
this  type  of  flow  is  significantly  more  complicated  due  to  the 
close  proximity  of  several  fixtures  which  absorb  or  reflect 
sound  in  varying  degrees.   In  a  turbojet  test  cell  the  engine 
is  tested  in  close  proximity  to  an  augmentor  tube  which  is 
used  for  mixing  and  cooling  the  exhaust  gases  before  they  are 
vented  to  the  atmosphere  through  a  stack.   Each  of  these  struc- 
tures influences  the  sound  propagation  in  some  way. 

The  numerical  models  discussed  above  can  be  used  to  predict 
velocity,  pressure  and  turbulence  intensity  distributions  in 
confined  turbulent  flows  like  those  discussed  above.   It  may 
be  possible  to  use  the  computer  codes  to  predict  optimum  place- 
ment of  sound  suppressors  in  geometrically  complicated  flows 
like  those  in  a  turbojet  test  cell. 

An  attempt  was  made  to  measure  the  sound  produced  within 
the  Naval  Postgraduate  School  subscale  turbojet  test  cell  aug- 
mentor tube.   These  data  were  compared  to  the  flow  field  pre- 
dictions obtained  with  the  computer  model. 
B.   EXPERIMENTAL  APPARATUS 

1.   Test  Cell 

Design  and  construction  of  the  NPS  subscale  turbojet 
test  cell  are  detailed  in  reference  15.   The  cell  is  a  one- 
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eighth   scale   model   of   a  NAS   Alameda   test  cell.      A   TF   41 
engine  was   scaled  to  one-eighth   in   diameter  resulting  in  the 
mass    flow  being   scaled  by  one-sixty    fourth    to   maintain    flow 
velocities   equal   to   those   in   the   full   scale   cell.      The   test 
section  is   enclosed  by  hinged  plexiglas     aides    to   allow  easy 
access   and  visual  monitoring  of  the  section  during  operation. 
These   sides  were  held  open  during  the   acoustic  measurements 
described   in    this    report.      The   augmentor   tube,    equipped  with 
a   straight  inlet,    exits    the   cell    through   a   removable  wall. 
Its    downstream  end  is   normally   attached   to   a   deflector-plate- 
equipped  vertical  exhaust  stack.      For  this   experiment,   however, 
the   stack  was    removed, allowing   the    fluid  exiting  the   augmentor 
tube    to   vent    to    the   atmosphere.      In    an  effort   to    reduce    sound 
reflection   inside    the    test   cell,    a   one-quarter   inch  plywood 
divider  with    three-quarter   inch    foam  padding  was   extended    from 
the   augmentor   inlet    to    the   plexiglas    sides   of   the   test   cell. 

The   combustor  used   to   simulate    turbojet   tailpipe    and  nozzle 
conditions   was    supplied  with    compressed   air    from  an  Allis- 
Chalmers,    twelve    stage    axial    compressor.      The    tailpipe  was 
fitted  with   either   a   small    (2    inches    diameter)    or   large    (3.5 
inches    diameter)    nozzle. 

2 .      Ins  trumentation 

Several   methods   were   employed   in   an   attempt   to   measure 
the   aerodynamic   noise   generated  by    the   engine   exhaust  in    the 
subscale    test   cell.      The    proximity   of   various    structures    in- 
cluding   the   augmentor   tube    and   the    stack   produced   an    "acoustic- 
ally   dirty"    environment.      This   prevented   useful    data    from  being 
taken   external    to    the    augmentor   tube    regardless   of    the   pressure 
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transducer  orientation  to,  or  distance  from  the  jet  axis.   The 
method  chosen  to  obtain  these  measurements  is  illustrated  in 
figure  2. 

The  top  of  the  augmentor  tube  incorporated  a  series  of 
pressure  taps.   Several  of  these  were  modified  to  accept  a 
one-eighth  inch  directional  microphone.   The  microphone  was 
inserted  flush  with  the  inside  of  the  augmentor  tube.   Ini- 
tially it  was  held  in  a  plexiglas   plug  with  a  set  screw 
(figure  2) .   The  transducer  signals  were  amplified,  displayed 
and  recorded  using  a  Spectral  Dynamics  Corporation  model  50  330 
spectrum  analyzer  and  a  HP  70  35B  X-Y  recorder. 

After  initial  testing  an  attempt  was  made  to  isolate  those 
frequencies  associated  with  the  physical  resonance  of  the  aug- 
mentor tube  and  test  cell  from  those  generated  directly  from 
the  jet  mixing  process.   The  pressure  transducer  was  wrapped 
with  styrofoam  and  refitted  in  the  plexiglas   holder  without 
the  securing  set  screw.   To  determine  the  degree  to  which  the 
adjacent  building  was  contributing  to  the  recorded  data,  a 
plywood  sheet  was  positioned  between  the  cell  and  building. 
C.   DATA  ACQUISITION  AND  REDUCTION 

Sufficient  data  were  recorded  for  each  test  condition  to 
allow  the  mass  flow  rates  for  the  nozzle,  test  cell  and  aug- 
mentor tube  to  be  calculated.   In  addition,  a  plot  of  frequency 
vs  sound  intensity  was  recorded  at  each  of  the  seven  augmentor 
ports.   On  subsequent  runs  with  the  modified  transducer  holder 
and  alternate  background,  data  were  recorded  under  similar 
conditions  at  augmentor  port  number  four  only. 

16 


D.   COMPUTER  MODEL  PREDICTIONS 

At  attempt  was  not  made  to  predict  the  sound  distribution 
within  the  augmentor  tube.   Rather  the  velocity,  pressure  and 
turbulence  kinetic  energy  distributions  were  calculated  using 
the  computer  program  discussed  above.   It  was  desired  to  deter- 
mine whether  these  distributions  could  be  used  as  a  first  ap- 
proximation for  locating  the  peak  noise  intensities  within  the 
augmentor  tube. 
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IV.   DISCUSSION  OF  RESULTS 

A.   INTRODUCTION  . 

The  purpose  of  this  study  was  to  utlize  a  primitive  varir- 
able,  finite  difference  computer  program  to  analyze  the  flow 
within  a  turbojet  test  cell  and  to  validate  the  model  with  data 
obtained  from  a  subscale  turbojet  test  cell  located  at  the 
Naval  Postgraduate  School.   Previous  work  [2,6,7]  had  accom- 
plished this  task  for  a  ty-ui     computer  model  and,  therefore, 
empirical  data  and  the  predictions  of  the   ip— oo  model  were  avail- 
able for  comparison.   The  experimental  data  available  consisted 
of  augmentor  wall  pressure  distributions  and  radial  velocity 
profiles  along  the  length  of  the  augmentor  tube  for  low,  medium, 
and  high  engine  flow  rates. 

As  previously  indicated,  the   i^-co   model  did  a  poor  job  of 
predicting  pressure  distributions  in  all  but  quiescent  flow  re- 
gions.  Additionally,  numerical  convergence  was  difficult  to  ob- 
tain with  that  model  when  used  to  predict  high  velocity  flows 
where  compressibility  effects  are  significant.   It  was  antici- 
pated that  a  primitive  variable  model  would  help  to  alleviate 
these  difficulties.   It  is  desirable  to  have  a  model  which  can 
be  used  to  predict  the  flow  field  for  full-throttle  engine  con- 
ditions where  the  engine  exhaust  flow  is  choked. 

Adaptation  of  the  primitive  variable  model  to  the  subscale 
test  cell  geometry  required  the  use  of  several  approximations: 

(1)   In  modeling  axi-symmetric  flow,  the  engine  was  by 
necessity  positioned  at  the  axis  of  symmetry.   In  the  actual 
test  cell  the  engine  was  mounted  closer  to  the  deck  than  to 
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the  overhead.   It  would  be  expected,  therefore,  that  the 
velocity  distribution  in  the  secondary  flow  (the  flow  around 
the  engine)  would  be  somewhat  different  than  predicted. 

(2)  The  subscale  test  cell  cross  section  was  rectangular 
while  the  engine  and  augmentor  tube  were  cylindrical.   The 
system  was  modeled  as  three  concentric  cylinders  with  cross- 
sectional  areas  equivalent  to  the  physical  system.   The  nozzle 
exit  area,  the  test  section  cross-sectional  area  and  the  empir- 
ical augmentation  ratios  and  mass  flow  rates  were  used  in  cal- 
culating the  axial  inlet  velocities  used  in  both  computer  models, 

(3)  The  actual  engine  incorporated  a  converging  nozzle. 
The  engine  was  modeled  as  a  cylinder  with  a  diameter  equal  to 
the  actual  nozzle  exhaust  diameter. 

(4)  In  the  model  the  augmentor  inlet  and  the  aft  test 
cell  wall  were  taken  to  be  flush.   The  actual  augmentor  is  often 
inserted  into  the  test  section  which  forms  a  recirculation  zone 
above  the  augmentor  tube.   The  effects  on  the  augmentor  flow 
field  introduced  by  this  recirculation  region  have  been  shown 

to  be  minimal  [7]  .   It  should  be  noted  that  the   i^-co  model 
incorporated  this  recirculation  zone  and  a  flow  reducing  lip 
flange  on  the  augmentor  inlet.   When  comparisons  were  made  bet- 
ween the  predictions  of  the  two  models,  the  effects  of  this  re- 
circulation zone  and  the  augmentor  inlet  lip  flange  in  the   ^-w 
model  were  minimized  by  reducing  their  dimensions  to  one  grid 
spacing. 

There  was  one  additional  difference  between  actual  test 
conditions  and  the  numerical  model.   In  a  test  cell  of  speci- 
fied geometry  (cross-sectional  area,  inlet  and  stack  acoustic 
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and  pollution  treatment  devices,  augmentor  tube  length,  dia- 
meter and  inlet  geometry,  etc.)  the  mass  flow  rate  of  secon- 
dary air  pumped  by  the  engine  exhaust  and  the  pressure  distri- 
bution throughout  the  cell  are  unique  functions  of  the  engine 
flow  rate.   The  static  pressures  at  the  cell  inlet  and  stack 
exhaust  are  both  equal  to  the  local  atmospheric- pressure. 
Thus,  the  most  desirable  model  capability  would  be  the  predic- 
tion of  the  cell  flow  field  as  a  function  only  of  engine  flow 
rate.   To  accomplish  this  with  a  two-dimensional  finite-differ- 
ence model  would  greatly  complicate  the  analysis  and  increase 
the  required  computer  computation  time.   In  this  investigation 
the  flow  field  was  predicted  only  for  the  test  cell  flow  from 
the  engine  inlet  to  the  exhaust  of  the  augmentor  tube.   Since 
inlet  boundary  conditions  are  required  for  the  solution  of  the 
elliptic  equations,  the  mass  flow  rate  of  secondary  air  was 
specified  at  the  experimentally  measured  valve.   The  predicted 
pressure  profiles  were  then  compared  to  the  measured  values. 
An  alternate  method  could  have  been  employed,  requiring  signi- 
ficantly greater  computer  time:   iterate  on  the  inlet  mass 
flow  rate  until  a  specified  augmentor  exhaust  pressure  was 
attained. 

Three  of  the  flow  conditions  selected  for  model  valida- 
tion corresponded  to  conditions  where  empirical  data  were 
available.   Model  predictions  for  two  additional  conditions 
were  made  increasing  the  engine- to-augmentor  inlet  spacing  to 
one  and  two  engine  diameters.   No  empirical  data  were  avail- 
able for  the  two  latter  conditions  and,  therefore,  inlet  para- 
meters were  simulated  using  empirical  data  for  zero  engine- 
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to-augmentor  entrance  spacing.   The  test  conditions  are 
summarized  in  Table  III. 

Figures  3  through  7  compare  predicted  axial  pressure 
distributions  and  radial  velocity  profiles  obtained  with  the 
ij/-u)  and  u-v-p   computer  models.   In  addition,  the  available 
empirical  data  are  also  plotted  on  those  figures.   The  velocity 
profiles  from  both  computer  models  were  plotted  for  the  grid 
lines  closest  to  the  locations  of  the  experimental  data.   In 
the  cases  where  empirical  data  were  not  available,  various 
representative  velocity  profiles  were  plotted  for  both  models. 
Experimental  pressure  profiles  were  available  only  on  the  deck 
of  the  test  cell  and  along  the  top  of  the  augmentor  tube.   Thus, 
for  negative  positions  in  Figures  3  through  7  no  empirical 
data  were  available  for  direct  comparison  to  the  predicted  pro- 
files.  The  models  did  predict  a  very  small  pressure  drop  along 
the  test  cell  deck  in  agreement  with  experiment.   Predicted 
axial  pressure  profiles  are  presented  for  various  radial  posi- 
tions.  These  locations  are  given  as  fractions  of  augmentor 

tube  radius   (R  )  .   For  example,  the  pressure  distribution 
a 

labeled   R  =  0.96  R   indicates  that  the  distribution  is  along 

a 

an  axial  grid  line  located  at  a  distance  96  percent  of  the  aug- 
mentor radius  from  the  axis  of  symmetry.   Two   ^-co  model  axial 
pressure  profiles  are  depicted  for  each  test  condition.   One 
profile   (R  =  0.38  R  )   lies  in  the  quiescent  flow  region  bet- 
ween  the  engine  and  augmentor  and  was  previously  found  by 
Walters  [6]  to  be  the  only  location  which  produced  a  reasonable 
estimate  of  the  measured  profiles.   The  second  axial  pressure 
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distribution   (R  =  0.13  R  )   was  located  at  less  than  one 

a. 

engine  radius  of  the  center  line.   Three   u-v-p  model  axial 

pressure  profiles  are  presented  for  each  test  condition.   One 

distribution   (R  =  0.04  R  )   is  near  the  axis  of  symmetry. 

a 

Another  profile   (R  =  0.28  R  )   runs  along  the  top  of  the 
engine  and  through  the  turbulent  mixing  region.   The  third 
profile   (R  =  0.96  R  )   is  close  to  the  augmentor  wall. 
B.   RESULTS  AND  DISCUSSION  -  PRESSURE  AND  VELOCITY  DISTRIBUTIONS 
1.   Test  Case  I  -  Low  Flow  Rate/Zero  Spacing 

a.  Velocity  Profiles 

Both  models  predicted  virtually  identical  velocity 
profiles  at  each  station  along  the  augmentor  tube  (Figure  3) . 
There  was  close  agreement  between  the  predicted  and  the  experi- 
mental velocity  profiles  at  the  exit  of  the  augmentor.   The 
latter  result  was  expected  since  the  velocity  profile  had  be- 
come fully  developed  near  the  augmentor  exit. 

b.  Pressure  Profiles 

(1)   i^-oj   Model,   This  test  condition  had  the 
lowest  flow  rates  and  offered  the  best  chance  for  agreement 
with  experiment.   The  initial  pressure  drop  and  the  pressure 
rise  were  underpredicted  in  the  outer  flow  region   (R  =  0.38  R  ) 
Nearer  the  axis  of  symmetry   (R  =  0.13  R  )   the  initial  pressure 
decrease  was  significantly  overpredicted,  but  the  magnitude 
(from  minimum  to  maximum)  and  profile  of  the  pressure  rise  were 
in  good  agreement  with  the  experimental  wall  profile.   The  pre- 
dicted pressure  curves  leveled  off  earlier  than  the  experimental 
data.   The  difficulty  in  obtaining  good  pressure  profiles  with 
this  model  are  evident.   Near  the  augmentor  exit  there  should 
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be  negligible  radial  pressure  variations.   The  solution  was 

converged  in  all  dependent  variables  (Table  V)  and  negligible 

pressure  variation  was  calculated  in  the  radial  direction  near 

the  augmentor  exit.   However,  the  pressure  profiles  along  the 

two  axial  paths  (calculated  by  summing   (3  p/8  x)  .  Ax   from  the 

inlet  end)  did  not  become  equal  near  the  augmentor  exit.   The 

latter  resulted  from  the  large  errors  in  the  predicted  profiles 

near  the  augmentor  inlet  (i.e.,   ^-§  in  the  expression 

oj  x 

p  -  p^.^^  =   /_,  (■*£)  Ax   is  very  inaccurate  near  the  augmentor 
x    rer      n      dx 
x=0 

inlet)  .   This  difficulty  with  the   i|;-oj   model  has  been  previous- 
ly discussed  by  Walters  and  Netzer  [6], 

(2)   u-v-p  Model.   Both  the  pressure  drop  and  rise 
within  the  augmentor  were  underpredicted.   However,  the  primi- 
tive variable  model  more  accurately  predicted  the  location  of 
the  minimum  pressure  and  seemed  to  level  off  at  approximately 
the  same  augmentor  position  as  the  experimental  data.   Walters 
and  Netzer  [6]  have  previously  shown  that  the  models  predict 
that  mixing  is  nearly  complete  at  the  location  where  the  pres- 
sure profile  levels  off  (and  simultaneously,  where  the  velocity 
distribution  approaches  a  fully  developed  profile) .   Experi- 
mental data  confirmed  this  characteristic.   Figure  3  indicates 
that,  as  anticipated,  the  u-v-p  model  can  more  accurately  pre- 
dict the  location  at  which  turbulent  mixing  is  complete.   It 
also  predicted  very  little  pressure  variation  with  radial 
augmentor  position  as  is  known  to  be  true  experimentally. 
2 .   Test  Case  IIA  -  Medium  Flow  Rate/Zero  Spacing 
a.   Velocity  Profiles 

The  results  are  presented  in  Figure  4  and,  again, 
both  computer  predictions  were  very  similar  and  agreed  with 
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the  limited  experimental  data  at  the  augmentor  exit  plane, 
b.   Pressure  Profiles 

(1)  i/j-u)   Model..   The  pressure  profile  nearest 
the  center line  became  more  unrealistic  for  this  higher  flow 
rate  condition.   The  initial  pressure  drop  was  greatly  exag- 
gerated; however,  the  magnitude  of  the  pressure  rise  was  again 
in  good  agreement  with  experiment.   The  pressure  profile  in 
the  quiescent  flow  region   (R  =  0.38  R  )   did  not  agree  with 

3. 

the  experimental  curve.   Both  the  pressure  drop  and  rise  in 
the  augmentor  tube  were  underpredicted.   In  addition,  the 
minimum  pressure  point  was  predicted  to  occur  about  one  engine 
diameter  downstream  of  the  experimental  minimum. 

(2)  u-v-p  Model.   The  u-v-p  model  appears  to 
consistently  underpredict  the  augmentor  pressure  drop  and 
rise.   The  slope  of  the  pressure  rise,  however,  was  in  reason- 
able agreement  with  experiment.   As  discussed  above,  the  pre- 
dicted pressure  profiles  were  determined  by  specifying  the 
inlet  mass  flow  rates  and  uniform  velocity  profiles  at  the 
engine  exhaust  and  cell  inlet.   Table  IV  presents  the  effects 
of  changes  in  the  specified  flow  rates  on  the  predicted  maxi- 
mum pressure  rise.   Also  shown  in  Table  IV  are  the  correspond- 
ing values  predicted  using  the  one-dimensional  equations  for 
conservation  of  mass,  momemtum  and  energy.   The  latter  calcu- 
lations were  made  assuming  perfect  gases,  coplanar  engine  exit 
and  augmentor  inlet  planes,  no  losses,  and  plug  flow  for  both 
inlets.   The  1-D  calculations  do  not  correspond  directly  to 
the  2-D  calculations  (no  losses  and  plug  flow  of  secondary  air 
at  augmentor  inlet) .   However,  the  predicted  magnitudes  of 
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Ap     and  the  effects  of  changing  flow  rates  should  be  similar, 
max 

Table  IV  indicates  that  the  2-D  and  1-D  model  predictions  for 
the  nominal  case  (II  A)  were  in  good  agreement.   However,  the 
experimental  APmax  was  considerably  larger  than  predicted. 
Cases  II  E  and  II  F  show  the  predicted  effects  of  10%  changes 
in  the  engine  and  secondary  flow  rates,  respectively.   The 
predicted  1-D  effects  were  less  than  the  2-D  effects  but  in  the 
correct  direction  and  sequence.   It  is  seen  that  a  net  10%  error 
in  the  measured  augmentation  ratio  (m   from  ASME  orifice, 
mv   from  augmentor  exit  pi  tot  rake,  and  m  =  m  -  m  )   could 
produce  the  differences  between  the  experimental  and  2-D  model 
predictions.   For  example,  an  m   larger  by  3%  and  a  total 
augmentor  flow  rate  smaller  by  3%  (which  are  within  the  experi- 
mental accurary)  would  produce  a  decrease  in  augmentation  ratio 
of  approximately  8%.   Thus,  the  differences  between  the  pre- 
dicted and  measured  pressure  profiles  may  not  be  as  large  as 
depicted  in  Figure  4. 

3.   Test  Case  IIB  -  Medium  Flow  Rate/One  Engine 
Diameter  Spacing 

a.   Velocity  Profiles 

Experimental  data  were  not  available  for  this  test 
condition.   The  predicted  profiles  are  presented  in  Figure  5 
and  show  that  the  two  computer  models  predicted  velocity  pro- 
files which  were  in  close  agreement.   The  primitive  variable 
model,  however,  predicted  a  slightly  greater  initial  jet 
spreading  at  the  augmentor  entrance  and  required  a  slightly 
longer  duct  length  to  obtain  a  fully  developed  profile. 


27 


b.   Pressure  Profiles 

(1)  ip-co   Model.   As  for  Case  IIA,  the  centermost 
pressure  distribution  greatly  exaggerated  the  initial  pressure 
drop.   In  this  case,  the  pressure  rise  also  appeared  to  be 
much  too  rapid  in  comparison  to  the   R  =  0.38  R   profile 

and  the   u-v-p  profiles.   The  profile  in  the  quiescent  flow 
region  was  in  reasonable  agreement  with  the  wall  profile  ob- 
tained with  the   u-v-p  model.   The  slopes  of  the  pressure 
rise  and  the  minimum  pressure  predicted  by  both  models  were 
nearly  identical.   The   ^-00   model,  however,  again  predicted 
the  minimum  pressure  to  occur  somewhat  farther  downstream  than 
did  the   u-v-p   model. 

(2)  u-v-p  Model.   For  this  case,  the  pressure 

profile  closest  to  the  augmentor  wall  had  a  significantly 

lower  minimum  pressure  than  the  other  profiles. 

4.   Test  Case  IIC  -  Medium  Flow  Rate/Two  Engine 
Diameter  Spacing 

a.  Velocity   Profiles 

No   empirical   data  were    available    for   this    test 
condition.      The    computer  predictions    are   presented   in   Figure    6. 
As    for   Case    IIB,    the   primitive   variable   model   predicted   greater 
initial    jet   spreading.       In   this    case,    however,    the      ip-oo 
velocity  profiles   became    flat   and   the   pressure   profiles    leveled 
off   considerably   upstream  of   the    u-v-p      predicted  profiles. 

b.  Pressure   Profiles 

(1)  ^-oo  Model-  The  centermost  pressure  profile 
was  completely  unrealistic.  The  quiescent  region  profile  in- 
dicated a  larger  pressure  drop  and  rise  and  a  minimum  pressure 
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point  farther  downstream  than  the  u-v-p  model.   In  addition, 
the  ip-co  model  profile  leveled  off  much  earlier. 

(2)   u-v-p  Model.   Again  the  minimum  pressure  was 
obtained  for  the  profile  closest  to  the  augmentor  wall  and 
was  located  just  inside  the  augmentor  entrance.   All  the  primi- 
tive variable  curves  indicate  a  more  gentle  augmentor  pressure 
rise,  leveling  about  midway  down  the  augmentor  tube. 
5.   Test  Case  III  -  High  Flow  Rate/Zero  Spacing 
a.   Velocity  Profiles 

Substantially  more  experimental  data  were  avail- 
able for  this  test  condition.   In  the  experiment  the  nozzle 

was  operated  with  a  pressure  ratio   (p  _/pm)   less  than  cri- 

atm  T 

tical  for  one-dimensional  isentropic  flow.   However,  using 
the  experimental  nozzle  flow  rate  and  approximating  the  nozzle 
flow  as  one-dimensional  and  isentropic  resulted  in  a  nozzle 
exit  Mach  number  of  approximately  0.95.   This  condition  was 
imposed  on  the  models.   As  has  been  observed  for  the  previous 
test  cases,  the  velocity  profiles  for  both  models  were  quite 
similar  (Figure  7)  and  in  reasonably  good  agreement  with 
experiment.   However,  the  predicted  profiles  for  both  models 
tended  to  flatten  a  little  too  rapidly.   The  primitive  vari- 
able model  again  predicted  slightly  more  initial  jet  spread- 
ing at  the  augmentor  inlet  and  less  mixing  down  the  augmentor 
tube  than  did  the  \p-u      model.   The  experimental  data  more 
nearly  agreed  with  the   i|>-oj   model  at  the  augmentor  inlet  and 
with  the  u-v-p  model  downstream. 
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b.   Pressure  Profiles 

(1)  iJj-u)  Model.   Again  the  center  pressure 
profile  was  completely  unrealistic.   The  quiescent  region  pres- 
sure profile  underpredicted  the  augmentor  pressure  drop  and 
rise.   The  minimum  pressure  position  was  substantially  dis- 
placed down  the  augmentor  tube  and  the  slope  of  the  pressure 
rise  did  not  agree  with  experiment. 

(2)  u-v-p  Model..   All  three  u-v-p  pressure 
profiles  were  quite  similar.   They  underpredicted  both  the 
augmentor  pressure  drop  and  rise.   The  slope  of  the  pressure 
rise  was  in  good  agreement  with  experimental  data  within  the 
first  half  of  the  augmentor  tube.   The  predicted  pressure  pro- 
files from  both  the   ip— co   and  u-v-p  models  leveled  off  before 
the  experimental  curves. 

C.   COMPUTATIONAL  ACCURACY  AND  REQUIRED  COMPUTER  TIME 

The  utility  of  any  computer  program  using  numerical 

methods  is  reflected  by  the  amount  of  CPU  time  required  and 

the  ease  of  arriving  at  converged  solutions.   Table  V  compares 

the  percentage  change  in  variable  magnitude  on  successive 

iterations  and  the  required  CPU  time. 

A  considerable  savings  in  CPU  time  was  obtained  using 

the  line-by-line  iterative  procedure  of  the  primitive  variable 

model  in  lieu  of  the  point-by-point  (Gauss-Seidel)  method 

employed  in  the   ip— co  model. 

At  low  flow  rates,  the  convergence  was  quite  similar  for 

both  models.   However,  at  higher  flow  rates  the  u-v-p  model 

had  better  convergence  in  less  time. 
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D.   COMPUTER  RELATED  PROBLEMS 

Both  models  required  that  the  proper  relaxation  para- 
meters be  selected  in  order  to  obtain  convergence.   The  lack 
of  any  procedure  for  selecting  the  proper  relaxation  values 
makes  this  process  quite  time  consuming.   Previous  research 
using  the   ip-co  model  facilitated  the  selection  of  these 
parameters  for  that  model.   It  was  found  that  the  u-v-p  model 
was  quite  sensitive  to  the  calculated  "pressure  corrections". 
Obtaining  the  correct  underrelaxation  value  for  pressure 
proved  to  be  the  key  in  arriving  at  a  converged  solution  for 
the  primitive  variable  model. 

The  line-by-line  iterative  procedure  used  in  the  u-v-p 
model  was,  as  expected,  quite  good  in  propagating  distur- 
bances downstream  when  sweeping  from  left  to  right.   A  down- 
stream disturbance  is  propagated  upstream  by  successive  sweeps 
through  the  entire  field.   This  fact,  at  least  for  a  geometry 
incorporating  a  sudden  contraction,  made  the  convergence  de- 
pendent on  the  number  of  traverses  on  each  radial  line.   An 
excessive   number  of  traverses  would  cause  divergence.   The 
number  of  traverses  on  each  line  was  controlled  in  two  ways. 
After  each  traverse,  residual  factors  were  calculated  for 
each  variable  and  the  largest  residual  factor  was  compared  to 
a  pre-set  value.   If  the  largest  residual  was  less  than  the 
pre-set  value,  the  program  advanced  to  the  next  radial  line. 
The  program  would  also  advance  when  a  pre-set  maximum  number 
of  traverses  had  been  completed  on  any  radial  line.   To  aid 
convergence,  only  a  few  traverses  on  each  line  were  allowed 
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until  several  sweeps  through  the  field  had  allowed  the  pre- 
sence of  the  contraction  wall  to  be  "felt"  upstream.   The 
number  of  traverses  on  a  line  was  then  increased.   It  was 
additionally  found  that  when  working  with  coaxial  flows  with 
radically  different  inlet  velocities,  the  normalizing  factors 
(which  were  based  on  average  inlet  conditions  and  used  to 
calculate  the  residual  values  on  each  line)  resulted  in  ex- 
cessive traverses  being  made  in  regions  of  high  flow  velocity. 
Repeated  calculations  on  radial  lines  which  had  already  con- 
verged often  caused  divergence.   Adjusting  the  normalizing  fac- 
tors downstream  of  the  engine  exit  alleviated  this  problem. 

The  primitive  variable  model  demonstrated  some  conver- 
gence difficulty  in  the  recirculation  region  adjacent  to  the 
sudden  contraction.   This  problem  could  have  resulted  from 
the  relatively  large  normalizing  factors  used  in  this  local 
region  of  low  velocity  or,  as  suggested  by  Launder  and 
Spalding  [12],  it  could  have  possibly  been  due  to  the  in- 
adequacy of  the  empirical  constants  (Table  I)  in  the   K-e 
model  for  the  flow  in  this  quiescent  zone. 

As  with  any  finite  difference  numerical  solution,  grid 
spacing  was  found  to  be  critical.   To  aid  convergence,  the 
grid  spacing  was  decreased  in  all  regions  where  large  pro- 
perty gradients  were  expected.   Gosman,  et  al  [1]  recommended 
that  for  the   ^-os  model,  successive  spacing  should  not  in- 
crease by  more  than  about  a  factor  of  1.5.   This  restriction 
was  also  employed  for  the  primitive  variable  model.   A  30  by 
30  grid  system  was  utilized  for  the  primitive  variable  model, 
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whereas  a  4  3  by  4  0  grid  was  required  for  the   i^-oo  model. 

E.   SUMMARY  OF  RESULTS  -  PRESSURE  AND  VELOCITY  DISTRIBUTIONS 

In  most  cases  the  predicted  velocity  profiles  for  both 
the   ip-to   and  u-v-p  computer  models  agreed  with  each  other  and 
with  the  available  experimental  data.   The  u-v-p  model  predic- 
ted more  initial  jet  spreading  at  the  augmentor  inlet,  espe- 
cially in  those  cases  where  the  engine  exit  plane  was  not  flush 
with  augmentor  inlet.   Both  models  predicted  that  a  flat  velo- 
city profile  was  obtained  where  the  pressure  profile  leveled 
off.   For  the  one  case  in  which  experimental  data  were  avail- 
able over  the  entire  legnth  of  the  augmentor  tube,  the  pre- 
dicted velocity  profiles  seemed  to  flatten  slightly  faster 
than  the  experimental  data. 

In  general,  the   i|>-oo  model  demonstrated  poor  pressure 
prediction  capability.   For  low  flow  rates  the  pressure  rise 
on  the  wall  from  minimum  to  maximum  was  accurately  predicted 
for  a  pressure  profile  calculated  near  the  axis  of  symmetry. 
For  higher  flow  rates  the  centermost  predictions  became  un- 
realistic although  converged  solutions  were  obtained  for  all 
primary  variables.   For  the  predicted  pressure  profiles  along 
axial  lines  that  were  in  the  quiescent  flow  region,  the  ip-w 
model  characteristically  unpredicted  the  pressure  drop  and 
rise  along  the  augmentor  wall.   The  minimum  predicted  pres- 
sure location  was  typically  somewhat  downstream  of  its  experi- 
mentally determined  position. 

The  primitive  variable  model  predicted  little  radial  pres- 
sure variations  except  near  the  engine  exit  plane.   It  seemed 


34 


to  accurately  locate  the  minimum  pressure  position  and  to  pre- 
dict the  rapid  pressure  drop  at  the  augmentor  entrance.   It 
consistently  underpredicted  the  pressure  drop  and  rise  in  the 
augmentor  but  tended  to  accurately  predict  the  slopes  of  the 
rising  pressure  profiles.   However,  as  discussed  above,  reason- 
ably expected  experimental  errors  could  account  for  these  dif- 
ferences.  Application  of  the  K-e   turbulence  model  with 
fixed  parameters  (C_,  C.  ,  C~,    a  ,  a„)   to  the  test  cell  flow 

U    X     Z  E    J\ 

conditions  may  also  be  a  reason  for  the  lack  of  quantitative 
accuracy  in  the  predicted  pressure  profiles.   Further  experi- 
mental data  (for  example,  turbulence  intensity  measurements) 
are  needed  to  further  validate  the  model. 

The  primitive  variable  model  appears  to  reasonably  pre- 
dict the  flow  field  up  to  sonic  engine  exhaust  conditions. 
These  predictions  include  realistic  pressure  distributions 
and  require  substantially  less  computer  time  and  fewer  grid 
lines  than  the   ip— oj  model. 
F.   RESULTS  AND  DISCUSSION  -  NOISE  CONSIDERATIONS 

Computer  model  predictions  for  the  velocity,  turbulence 
kinetic  energy  and  pressure  distributions  are  presented  in 
Figures  8,  9  and  10.   The  predictions  were  for  the  two 
inches  diameter  nozzle  at  high  mass  flow  rate  (Run  1,  Table 
VI) .   The  four  curves  on  each  figure  are  for  different  radial 
positions;  jet  centerline,  just  inside  the  jet  exhaust  radius; 
just  outside  the  jet  exhaust  radius,  and  near  the  augmentor 
wall. 
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Noise  intensities  for  selected  frequency  ranges  as  a 
function  of  augmentor  tube  position  are  displayed  in  Figures 
11,  12  and  13. 

Figures  14  through  16  and  Table  VII  compare  the  basic 
test  results  to  data  obtained  at  one  transducer  port  after 
the  microphone  was  insulated  and  the  "building  wall"  location 
was  changed. 

Figure  11  presents  the  data  for  the  small  nozzle  with  a 
high  subsonic  exhaust  velocity   (M  %  0.75).   The  highest 
frequencies  were  greater  than  peak  values  expected  from  ori- 
fice data   (V/2d)   but  were  apparently  generated  near  the 
nozzle  exit  as  expected.   If  the  high  frequencies  generated 
near  the  nozzle  exit  were  radiated  at  45°  to  the  axis,  the 
corresponding  peak  intensities  would  occur  between  two  and 
five  jet  diameters  downstream.   The  lowest  frequencies  peaked 
in  intensity  at  approximately  20  jet  diameters  downstream  in 
general  agreement  with  data  from  free  jets  emitted  from  orifice 
plates  [14].   The  intensity  increased  as  the  frequency  de- 
creased. 

Figure  12  presents  the  data  for  the  small  nozzle  with  a 
very  low  exit  velocity   (M   «  0.25) .   The  data  had  similar 
characteristics  to  that  for  the  high  exit  velocity:   (1)  the 
intensity  increased  with  decreasing  frequency,  (2)  the  peak 
frequencies  were  greater  than   V/2d  and  (3)  the  lowest  fre- 
quencies peaked  in  intensity  at  approximately  2  0  jet  diameters 
downstream.   In  this  case,  however,  the  highest  frequencies 
peaked  in  intensity  further  downstream. 
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Figure  13  presents  data  for  the  larger  nozzle  with  low 
exit  velocity.   The  data  were  similar  to  Figure  12  except 
that  the  intensity  did  not  vary  as  consistently  with  frequency 
range.   The  peak  intensity  did  occur  for  the  lowest  frequency 
range  and  occurred  at  approximately  10  jet  diameters  down- 
stream.  This  peak  intensity  was  at  the  same  augmentor  posi- 
tion as  for  the  smaller  nozzle  at  low  exit  velocity  (Figure  12) 

The  possibility  existed  that  some  of  the  sound  recorded 
was  either  being  reflected  from  or  was  due  to  the  physical 
resonance  of  surrounding  structures.   However,  the  data  pre- 
sented in  Figures  14  and  through  16  indicate  that  these  effects 
were  not  large. 

The  computer  predictions  made  for  the  distributions  of 
axial  velocity,  turbulence  kinetic  energy  and  pressure  (Figures 
8,  9  and  10)  were  for  the  test  conditions/results  presented  in 
Figure  11.   (RUN  No.  1,  Table  VI.) 

The  peak  intensities  for  the  lower  frequency  sound  noted 
above  occurred  approximately  at  the  predicted  location  within 
the  augmentor  tube  at  which,  (1)  turbulent  mixing  became  essen- 
tially complete,  (2)  pressure  reached  a  maximum,  and  (3)  the 
axial  velocity  profile  became  fairly  constant.   The  decay  in 
the  predicted  turbulence  kinetic  energy  (Figure  9)  followed 
closely  the  rise  in  noise  amplitude  for  the  two  lower  fre- 
quency ranges  (Figure  11) . 

The  data  were  much  too  limited  to  draw  any  general  con- 
clusions.  However,  the  results  indicate  that  the  maximum 
sound  intensities  occur  within  the  augmentor  tube  in  the 


39 


region  where  jet  mixing  is  nearing  completion.   If  this  ob- 
servation is  correct,  it  would  enable  the  computer  model  to 
be  used  for  full  scale  test  cells  to  predict  the  best  place- 
ment of  sound  suppressors/jet  breakup  devices   for  varying 
test  cell  geometries  and  engine  operating  conditions. 
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V.   CONCLUSIONS 
The  primitive  variable  computer  model  for  the  turbojet 
test  cell  appears  to  reasonably  predict  the  flow  field  for 
engine  exhaust  conditions  up  to  the  sonic  condition.   The 
conclusion  is  based  upon  subscale  test  cell  measurements  of 
velocity  and  pressure  distributions.   To  further  validate  the 
model,  comparisons  should  be  made  with  full  scale  test  cell 
data.   In  addition,  turbulence  intensity  distributions  mea- 
sured within  the  subscale  test  cell  would  provide  valuable 
new  data  for  additional  model  validation. 
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